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Abstract 



Averaged over ensenible of initial conditions kinetic transport equations of weakly coupled 
fH I systems of quarks and gluons are derived. These equations account for the correlators of fluc- 

tuations of particles and classical gluon fields. The isotropization of particle momenta by field 
fluctuations at the early prethermal stage of matter evolution in ultrarelativistic heavy ion colli- 
^ . sions is discussed. Our results can be useful for understanding under what conditions isotropiza- 

I tion of the quark-gluon plasma in ultrarelativistic heavy ion collisions can be reached within 

phenomenologically observed time scales. 
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^ . 1 Introduction 
00 . 

, Ideal fluid hydrodynamic models provide a good description of Relativistic Heavy Ion Collider 
(RHIC) data on single-particle hadron momentum spectra and elliptic flows [Ij. It is noteworthy 
that the agreement with data is achieved only if very rapid equilibration proper times Tj < 1 
fm/c are assumed f2|, whereas the theoretical estimates based on the initial conditions obtained 
I from the Color Glass Condensate (CGC) wave functions of colliding nuclei (see, e.g., Ref. [3j) and 
perturbative scattering processes yield thermalization times near 3 fm/c [4]. The possible resolution 
of this puzzle was proposed in Ref. [5] where it was argued that for hydrodynamic modeling of the 
early stage of RHIC collisions thermalization is not required and isotropization of parton momenta 
in the local fluid rest frame (local rest frame of the energy flow) suffices for applicability of ideal 
fluid hydrodynamic models. Then "early thermalization" in heavy ion collisions is more properly 
interpreted as evidence of the local isotropization in momentum space of the nonequilibrium quark- 
gluon plasma (QGP), which precedes the thermalization by acting on a faster time scale. The 
gauge fields' instability effects caused by a particle momentum anisotropy in the local rest frame (it 
appears due to the rapid longitudinal expansion), analogous to Weibel instability in Abelian plasmas 
[6], would speed up the onset of isotropization and subsequent thermalization of weakly coupled 
QGP (for review see, e.g., Refs. [Tl'H]) in relativistic heavy ion collisions because the unstable modes 
tend to make the particle momentum distributions more isotropic. Then a important question is 
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whether such instabihties could restore isotropy in momentum space in heavy ion colhsions on a 
relevant time scale. 

Recently the 3 + 1 dimensional numerical simulations of nonAbelian plasma instabilities in a 
stationary anisotropic weakly coupled plasma (g <^ 1 where g is the coupling constant) were carried 
out P [ini E] using the framework of the hard-loop (HL) effective theory jlS [l3l [H] . The HL 
effective theory treats the plasma particles as having arbitrary large momentum; this approximation 
is justified because in weak gauge coupling scenarios relevant for the early stage of relativistic heavy 
ion collisions |k| <^ |p|, where |k| is the characteristic momentum scale associated with the unstable 
gauge fields' modes and |p| is the typical momentum of plasma particles. This approach does not 
account for back reaction of the gauge fields on the particles because the gauge fields' effect on 
the particle "trajectories" becomes important when the amplitude of the gauge fields is ^ ~ IpI/s 
(in the Coulomb gauge), while in these simulations A < |k|/(7. Even with this simplification, the 
HL dynamics for nonAbelian theories are rather complex as a result of the nonlinear gauge field 
self-interactions that come into play when A ~ |k|/g. The simulations performed for a plasma 
with moderate momentum anisotropy [9] showed that unlike Abelian plasmas (and in contrast 
with earlier results of 1 + 1 dimensional simulations [15]) the gauge field dynamics changes from 
exponential field growth with time to linear growth when the vector potential amplitude reaches the 
nonAbelian scale, A ~ |k|/(5r, where nonlinear gauge field interactions become important. Evidently, 
it could slow down instability-driven particle momentum isotropization in heavy ion collisions. The 
same is valid for 3 + 1 dimensional numerical simulations [10] with extreme momentum anisotropy 
and strong initial fields that are nonperturbatively large, A > |k|/(7. As for the case of very 
strong momentum anisotropy and perturbatively weak initial fields, A <C |k|/(7, it was found [llj 
that the exponential growth continues beyond the nonAbelian bound and extends to higher wave 
vectors as compared to the perturbative scenario. The HL effective theory was also extended to 
the case of boost-invariant longitudinally expanding distribution of plasma particles [T6] and it was 
demonstrated in 1 + 1 dimensional simulations ^17j that chromo-Weibel instabilities grow nearly 
exponentially in the square root of proper time |18j . 

Important results have been also obtained beyond the HL techniques. The nonAbelian collective 
instabilities were studied in 3 + 1 dimensional numerical simulations [19] of Yang- Mills equations 
for unstable matter expanding into the vacuum after a high energy heavy ion collision. These 
calculations account for the back reaction of the soft field modes on the hard modes (" particles" ) . 
The calculations performed for the CGC initial conditions with fluctuations of the fields in rapidity 
(violations of boost invariance) demonstrated that nonAbelian self-interactions cause the growth of 
soft modes to saturate; however, the isotropization time scales of hard modes are much shorter if 
there are large initial fluctuations [19] . The solution of the full three-dimensional classical Vlasov 
transport equations [20] also goes beyond the HL approximation. In this kinetic approach, the mass 
shell partons released from the wave functions of relativistic colliding nuclei are treated as " fields" , 
if their momenta are much below the saturation momentum Qs (given by the square root of the 
color charge density per unit area in the incoming nuclei), and as "particles", if their momenta are 
on the order of Qs and above. It was found in numerical simulations of classical transport equations 
that for fairly strong initial random fields a very rapid isotropization of the particle momentum 
distribution was reached while there was no developed instability with rapid growth of the fields 
|20| . These results would indicate that for large initial field fiuctuations isotropization of particle 
momenta can be reached before the instabilities will develop. 

These findings probably indicate that thermalization of QGP at RHIC or Large Hadron Collider 
(LHC) energies is mainly an initial state problem and that fast isotropization in ultrarelativistic 
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heavy ion collisions can be reached if evolution starts from specific initial conditions. In this article, 
we demonstrate using analytical methods how fluctuations in the initial stage can speed up the 
isotropization process. For this purpose we, based on collisionless Vlasov transport equations of 
QGP [m |22l ia [M] and using the methods developed for Abehan plasmas [251 |26l El [28] (see 
also Ref. p9] where effective transport equations for nonAbelian plasmas are derived based on 
Wong equations [3^), derive kinetic equations of weakly coupled QGP (wQGP) that describe 
fluctuations driven isotropization of the averaged over ensemble of initial conditions distribution 
functions of on-mass-shell "particles" (quarks, antiquarks and gluons). The turbulent nonAbelian 
plasma instabilities are accounted for in an approach that was developed early for the turbulent 
Coulomb plasmas [3lj (for review, see Ref. [28j). 

This article is organized as follows. In Sec. 2, we derive averaged over ensemble of initial 
conditions kinetic transport equations of wQGP based on a collisionless approximation of kinetic 
equations. The "collision contributions", which naturally appear in these equations, are expressed 
through the correlators of fluctuations. In Sec. 3, we consider averaged wQGP kinetic equations 
below the nonAbelian scale where the nonlinear gauge fields' self-interactions can be ignored and 
we derive the corresponding collision terms that describe diffusion in momentum space leading to 
the locally momentum isotropical state as well as generalized Balesku-Lenard "collisions" leading 
to the (local) equilibrium state of the particle phase-space distributions. We conclude in Sec. 4. 

2 Averaged over ensemble of initial conditions kinetic transport 
equations of wQGP 

Let us start this section with a brief review of the transport equations of QGP in the collisionless 
Vlasov approximation. The collisionless approximation is applied to kinetic equations that describe 
the evolution of the distribution functions that are smoothed over physically small volume^ [25] , 
and this means neglect of short-range fluctuations that are responsible for the appearance of the 
Boltzmann collision terms in kinetic equations |25l [26] and, so, neglect of large angle ("hard") 
scattering of particles due to high transferred momenta at short distances. This approximation is 
justified for time scales that are short compared to the mean time between large angle scatterings 
of plasma particles. 

The distribution functions of quarks, Q{x,p), antiquarks, Q{x,p), and gluons, G{x,p), are 
assumed to satisfy the following collisionless nonAbelian Vlasov-type transport equations that 
describe high momentum modes that are treated as classical colored particles and soft gluons that 

^In a certain sense, utilization of "smoothed" quantities, and, so, transition from deterministic to probabilistic 
description, is unavoidable for macroscopic systems because it is impossible to fix ("observe") the micro-state of the 
macro-system with absolute accuracy without destroying-out the macro-state [32]. Note that even for exact initial 
conditions the "smoothing" appears effectively for the quantities that are calculated (in each "event") by means of 
numerical molecular dynamic models that are the "solver" of the reversible Hamiltonian dynamic equations. This is 
the result of the stochastic errors of the "round-up" of the numbers and the systematic errors of the method that are 
"triggered" by the dynamic chaos that is inherent to complex Hamiltonian systems. 
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are treated as classical fields (for details see, e.g., Refs. [71 \T3 \ \23 \ [Ml [33 | [M]): 

p^D,Q{x,p) + \p^{F,,, = 0' (1) 

p^^D, Q {x,p) - \p^{F,u, = 0' (2) 

p^D,G{x,p) + \p^{F,,, = 0- (3) 

Here x = {t,r), p = (po,p), 

D^ = Id^-ig[A^{x),...], (4) 

D^ = Td^-ig[\{x),...], (5) 

with and being four-potentials, F^jy and F^i/ being field strength tensors in the fundamental 
and adjoint representations, respectively, 

A, = Apa, (6) 

— ^'fi'^a, (7) 

Ff,^ = d^A^ - d^Af, - ig[A^, A^], (8) 
%u = d^A^- d^A^ - ig [A^ , 1^] , (9) 

where ta and are SU{Nc) group generators in the fundamental {Nc x Nc matrices) and adjoint 
[{N^ — 1) X i^c ~ 1) matrices] representations, respectively, a = 1, [N^ — 1), and the Einstein 
summation convention for repeated indices is utilized. The [...,...] denotes the commutator, and 
{...,...} denotes the anticommutator. 

Equations ([I|)-(l3]) are supplemented by the Yang-Mills equation 

D.F'^'^ix) = -fix), (10) 

where the color current density is expressed in the fundamental representation, = j^t", as 

f = ^[ d'p^{Qix,p)-Qix,p) - ^Tr[Qix,p) -Qix,p)]+2taTr[TaGix,p)]). (11) 

2 J Pq Nc 

Then, to derive kinetic equations for mean (statistically averaged) values one needs to perform 
ensemble average, (...), of Eqs. ([1]-([3|) and (jlOp over ensemble of initial conditions. The ensemble 
average allows split quark, antiquark, and gluon phase-space densities and gluonic classical fields 
into their mean part and a fluctuating part, for example, 

A^' = {A^") + 5A^'. (12) 

The mean value of the statistical fluctuations vanishes by definition, {6A^^) = 0. 

The quark and antiquark phase-space densities Q, Q are Nc x Nc hermitian matrices in color 
space and have singlet and multiplet parts in the fundamental representation of the SU (Nc) gauge 
group, and the gluon phase-space density G is an iN^ — 1) x (A'^ — 1) hermitian matrix in color 
space with singlet and multiplet parts in the adjoint representation. We assume that the statistically 
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averaged value of the partonic phase-space density is the singlet (and, so, is locally colorless) and 
the disturbance (fluctuation) of the phase-space density is the multiplet of the SU (Nc) gauge group. 
Then 

Q={Q) + 6Q = If, + 6f^ta, (13) 

Q=(Q)+SQ = If^+6f^ta, (14) 

G={G) + 6G = Ifg + 6f^Ta. (15) 

We assume that averaged values of the quark and antiquark phase-space densities coincide, fq = fg, 
and that the statistically averaged local value of the classical gluon field, A'^, and the gluon field 
strength, F^j,, are equal to zero, = (F^j,) = 0, and then ^4^ = and 

F^. = = ^M^^^ - 5-<5^M + 9SA';^SAtr''-. (16) 
Here we take into account that 



[t", t"] = i^H", [^^ T"] = z/"'«=r^ (17) 

where f""^^ are the antisymmetric SU{Nc) structure constants. Keeping quadratic terms in fluctu- 
ations, g6A'^6A'^f'^'^°', is necessary because one cannot ignore the nonAbelian gauge fields' self- 
interactions when corresponding amplitudes of fluctuations reach the nonAbelian scale where 
{d^SAf, — d^SAIj) ~ gdA'^^SA'^f^'^"'. Note here that the statistical average of the gluon field strength 
(stress tensor), {F^^,), is not only given by F^i,{{A)) due to quadratic terms in the fluctuations 
contained in F^y] therefore condition {Ffj_u) = implies 

{5A15AI) ~ 5'^\ (18) 

and, so, we conclude that the vanishing averaged local value of the classical gluon field strength 
means that fluctuations of different color components of classical gluon fields are statistically inde- 
pendent. 

Performing the ensemble average of Eqs. ([I|)-(l3]) and (jlOp gives 

p^{D^Q) + ^-p^{{5F^., ^}) = 0, (19) 



P'{D,Q) - \p^{{5F,,, ^}) = 0, (20) 
p^^D^G) + h^m,,, f^}) = 0, (21) 



and 



{D^dF^"") = 0. (22) 

For further convenience let us introduce mass-shell distribution functions fi{x,p) and 5fi{x,p) 
(index i here means q, q, or g) depending on four-position x = (t,r) and three-momentum p, e.g.. 



fq{x,p) = 2poe{po)d{p -m )fq{x,p) = 5{po-Jp'^ + ml)fq{x,p). (23) 
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For the adopted normalization 

j (fpfq{x,p)= j d^pfq{x,p) =nqix), (24) 

where nq[x) is quark number density. 

The distribution functions of quarks, antiquarks, and gluons have no simple probabilistic inter- 
pretation due to the gauge dependence since a color of a particle (e.g., a quark) can be changed 
by means of a gauge transformation. Only the traces of the distribution functions are gauge in- 
dependent and therefore they have the probabilistic interpretation (see, e.g., Ref. [33j). To get 
equations that govern the evolution of the colorless mean values and, so, present the average (most 
probable) evolution of the ensemble of systems, one can calculate Trace of Eqs. ((l9])-(l2l])El Be- 
cause Trace of the color multiplet part is equal to zero, the Trace explicitly reveals the evolution 
of the colorless (singlet) quantities, and the statistically averaged collision term does not change 
the color neutrality. Performing integration in Eqs. (jl9|) -()2ip over pQ and taking into account that 
Tr{t°') = Tr{T"') = and that Tr{AB) = Tr{BA) for any quadratic matrices of the same order, 
we gelH 

p^a,/,(x,p) = -j^p^Tr{{5F,^{xf-^^^)), (25) 
p'^dM^^v) = j^fTr{{5F,,{x)^-^^^)), (26) 
P^d.ax^p) = -Jf^fTr{{5F,^{xf-^^^)). (27) 

Equations (j25p - (l27p describe the evolution of the averaged over ensemble of initial conditions phase- 
space densities of quarks, antiquarks and gluons with collision terms that are determined by the 
correlators (statistically averaged products) of fluctuations. These "collision contributions" contain 
entire physics of "soft" (long distance) processes, e.g., isotropization in momentum space due to 
the chromodynamic Weibel instabilities and establishment of (local) equilibrium (thermalization) . 
Further analysis is limited to the quark collision term 

/,(x,p) ^ -j-p^Tr{{5F,^{xf-^^^)), (28) 

but expansion of the results to antiquark, Iq, and gluon, Ig, collision terms is straightforward. The 
detailed calculation of Iq is carried out in the next section for amplitudes of fluctuations below 
the nonAbelian scale; here we only demonstrate why the collision term, Iq, can lead to momentum 
space diffusion of particle phase-space density, /g, and, so, to local isotropization in momentum 
space. To demonstrate it, let us rewrite Iq in a slightly different form, 



Iq{x,p) = -j-±-Tr{{p^5F,,{x)5Q{x,p))) + Arr((5F^^(x)^^^^)). (29) 



Note here that the second term on the right-hand side of the above expression disappears in the 
Abelian approximation when, in particular, the stress tensor 5F^y is approximated by d^j^SAy — 



^The Trace of Eq. (I22J leads to trivial identity = 
^Starti 



^Starting from here the derivatives over po anywhere throughout the article are identically zero and po = Ep 
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dySA^. Then, subtracting Eq. (jl9p from Eq. ([T]), performing integration over pQ, and accounting 
for the leading order in fluctuations we get 

p^^D,5Q{x,v) = -gp^5F^,{x)^^^^ -igp^^{[5A,{x),5Q]). (30) 

OPu 

Here we take into account that [6A^{x),{Q)] = and that (1/2) ^} = SF^,y^. The 



formal solution of Eq. (j30p has the following form: 



SQix,p) = -gj d^yGpix - y)U{x,y)pnF^,{y)U{y,x)^^^^^^ + {other terms). (31) 

Here 

p^d^Gp(x)=5'^^\x), (32) 

Gp{x) = E~^Q{t)5^^\v - {p/Ep)t), (33) 

and U{x,y) is the gauge parallel transporter (see, e.g., Refs. [13\ 114 1 [23]) along the straight line 7 
going from y to x, 



U{x, y)=Vexp\-ig [ dz^ A^'{z)\ , (34) 
J7 



where V denotes path ordering. Note that in the Abelian approximation, when the terms that are 
not of the leading order in g are neglected, the transporter U{x,y) is approximated by unity. Then, 
substituting Eq. ([3T|) for 5Q{x,p) in the first term of Eq. ([29]) we get 

Ig{x,p) = j^-^ J d^Gpix - y)T,yf3{x,y)^^^^^^ + {other terms), (35) 

where 

%fs{x,y) = Tr{{p''5F^,{x)U{x,y)p^5F^p{y)U{y,x))). (36) 

Now, if we assume that the correlation lengths in {x — y) of T^p are far less than the time and 
length scales on which distribution function (Q) varies, we get 

Iq{x,p) = - — Dyp{x,p) — h {other terms), (37) 

OPu dpf3 

where 

F>uf3{x,p) = J d^yGp{x - y)%i3{x,y), (38) 

and we take into account that (Q) = fgl. Thereby, Ig contains a term that describes diffusion in 
momentum space of particle phase-space density and, therefore, Iq (as well as Ig, Ig) can be related 
with isotropization processes in relativistic heavy ion collisions. 



*We keep here non-linear terms in the covariant derivative and in the gluon field stress tensor because the corre- 
sponding nonAbelian terms cannot be neglected if the amplitudes of fluctuations reach the nonAbelian scale. 
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3 Kinetics of wQGP below the nonAbelian scale 

The main difficulty in calculating the collision terms is caused by the gauge fields' self-interactions 
that take place because of the nonAbelian nature of the QCD. In this section we discuss simple 
and theoretically clean situations when amplitudes of fluctuations are below the nonAbelian scale. 
Then one can neglect self-interactions and substitute Ffj_u — > = df^SA,^ — d^SA^^, Z)^ — > 9^. 

To obtain the explicit expressions for the correlators of fluctuations one needs to derive evolu- 
tional equations for the fluctuations. Then, subtracting averaged Eqs. ()19|) - (I22| ) from Eqs. ([I])-(l3l) 
and (jlOp . performing integration over pQ, and neglecting the terms that are not of the leading order 
in g, we get 

P>^d^5Q{x,p) = -5p/'5^^,(x)^Mi^, (39) 

p^d,6Q{x, p) = gpf^5T,Ax)^^^^, (40) 
p'^d^5G{x,p) = -gpf^6^^,(x)^^^^, (41) 



and 



d^dJ^^^'^ix) = -6fix). (42) 





= 5ff'^ 






J q 






= 5ff'^ 





The color current fluctuation, Sj'^{x), is expressed through fluctuations of the quark, antiquark and 
gluon phase-space densities and in the fundamental representation reads 

Sj^{x)=g J d'p^{6Q{x,p) - 5Q{x,p) +2taTr[TJG{x,p)]) = 

g J d'p^U^f^ix,p) - 6f^{x,p) + 2N,5f^{x,p)). (43) 

The solution of Eqs. (|39 p - (|4ip has the form 

(44) 
(45) 
(46) 

where Sfq^^\ 6f^^^\ and 5fg'^^^ are associated with the spontaneous fluctuations of the freely moving 
partons 

p''d^5f<^\x,p) =p^^d,5f^^'\x,p)=p^d,,6f<''\x,p) = 0, (47) 
and 6fq^^\ Sf^^^\ and 6fg^^^ are the partonic fluctuations induced by the classical field fluctuations 

SffHx, P) = -9 j d'^yGpix - y)p'^6J^-^{y)^Is^, (48) 

Sff\x,p)=g I d^yG^{x-y)p^5T;^{yf-^^l^, (49) 

5f<^{x,p) = -g j d^yG,ix-y)p>^6:F-^,iy)^^^. (50) 



8 



Then color current fluctuation can be presented as the sum of induced current fluctuation, (5j^(*) (x), 
and spontaneous current fluctuation, 6j^^^\x): 

mx) = Sj;:^'Hx)+6j!:^^\x), (51) 

where 

SjH^'Hx) = g J d3p^(5/^"(^)(x,p) - 5ff\x,p) + 2N,5f^^^\x,p)), (52) 

Sjf'> (x) = -g' J d^P^^d'yG,{x - y)p-5T^^{y)^i^, (53) 

and 

f{y, P) ^ fg{y, P) + MV, P) + 2NJg{y, p). (54) 



Substituting Eqs. (j5T]) - (j53|) into the right-hand side of Eq. (j42|) and moving the induced current 
fluctuation into the left-hand side, we get 

d^SJ^-^^ix) + 6r,^^\x) = -5j:^'\x). (55) 



Equation (|55|) can be solved in terms of four-potential fluctuations, 5A'^. The field SA'j^ is then the 

sum of the collective field fluctuations, dA'^j}^^"^''^ , that are the general solution of the homogeneous 
part of Eq. ([55]) . 



;;ions SA^ 

particle fluctuations. 



and the field fluctuations Jj4^*-^"'^*^ that are a particular solution of Eq. ([55|) and are related with 



Then, because 



we get 



^M5^a(part)(^) + 5j^;(^)(P'^^*)(x) = -5j<'\x). (57) 

5A1 (x) = 5A<^"^^'^ (x) + dAI^P'''-^'^ (x) , (58) 



5ff\x,p) = )(™")(x,p) + 5/,"W(f-*)(x,p), (59) 
5jf) (x) = <5j;^(*)M0 (x) + (x). (60) 

These formulas can be used to calculate correlators of fluctuations ("scattering terms") of 
averaged over initial conditions kinetic transport equations of wQGP in the Abelian approximation. 
Taking into account traceless of the stress tensor, 5J^li"'{x) = 0, one can rewrite the right-hand sides 
of Eqs. (j25p - (|27p in a slightly different form. 



p^a^/,(x,p) = -^^±.{pW^^(x)6f^{x,p)), (61) 
p>^d,Mx,p) = ^^±-y5r';,,{x)5f^{x,p)), (62) 



p^a^/,(x,p) = -^^^^±-y8T-^,{x)5rg{x,p)). (63) 
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Here we take into account that t", matrices satisfy 



TrlT^T^] 



1 



ab 



We present here details of the calculation of the scattering term, Iq 

g d 



[Ah) 



2N^ dp. 



{p^5K,{x)5f^{x,p)), 



(64) 
(65) 

(66) 



which governs the evolution of the statistically averaged quark distribution function in the Abelian 
approximation. The calculation of the other scattering terms is similar. 
Because of decomposition, Eqs. (HH) and ([59]) . one can see that 



(67) 



(5^»(f-*)(x)5/,"«(^'-*)(x,p)) + (5.f;(p-*)(x)5/,"W(x,p)). 



Let us assume that statistically averaged distribution functions are slowly varying in space and time, 
i.e., assume that the rate of evolution of the fluctuations is large compared to the inverse of time 
and length scales on which fq, fq, fg vary. To explicitly stress this assumption, we write space-time 
coordinates as index, e.g., fq(x){p) instead of fq{x,p) where necessary. Then, defining the Fourier 
transformations a{k) = f d^x ex.p{ikx)a{x) and a{x) = J ex.p{—ikx)a{k) for a function a{x), 
we get from Eqs. ([39])- dH]) for the k representation of the induced part of the particle phase-space 
distribution functions 



Sf:l^ik,p)=g 



6f;l^^{k,p) = -9 



<g!(fe,p) = <7 



-9^" + 



9'^" + 



9'"' + 



k^pf^ 


dfq{x)ip) 


Pak"" + iO 






dfq{x){p) 


Pak'' + iO 






dfg{x){p) 


Pak'' + iO 





-5Al{k), 
-SA^ik), 



and the k representation of the induced color current fluctuation 5j^^^^ is 

k^p" 



3 Of{x){p) 

^2E, dp. 



p^k"" + iO 



6A,{k) 



ng')(A:)M,(A:). 



Using Eq. (j7ip , we get the following expressiorj^l for Fourier transformed Eq. (j55p : 
(feV" - k^k'' - nf;)(A:))M,(A.) ^ e>l^^{k)6A,{k) = -5f^'\k), 

where the spontaneous color current fluctuation in fundamental representation, (k) = 6j^^^^ {k)ta 



(68) 
(69) 
(70) 

(71) 
(72) 



is Fourier transformed Eq. ([52^ : 
5j^'\k)=g [ d'p 



2E, 



■(5/«W(A:,p) - Sff\k,p) + 2NJf^('Hk,p))- 



(73) 



p 



^Note that the obtained equation is the same as that for Abelian plasmas, and if the right-hand side is assumed 
to be equal to zero the equation coincides with the linear response method equation for electromagnetic plasmas and 



wQGP (see, e.g., Refs. I71[T31[21IM1E1])- 



10 



Then the formal solution of Eq. (j72p as for four-potential fluctuations, 5A'^{k), is 



ef;)(fc)M«(-")(A;)=0, 



(74) 
(75) 

(76) 



Here the matrix (e ^)ufj.{x){k) is the inverse of ej^^(2. 

Now we can write the final expressions for fluctuations of the phase-space densities of quarks 
below the nonAbelian scale. Taking into account Eqs. ([H]) . ([7i|) . and ([Ml), we get 



Sf^^'\k,p)+g 



IIV 



Pak°' + iO 



-^1^^ ( A;) + M;^^^"'^*) ( A:) ) . 



(77) 



Then, using Eq. (j76j) . we get 



Now, using the relation 



9^" + 



9^" + 



Pak'^ + iO 



p'^,5.F^,(A;) = i-ip'^k^gi^ + ik^p^SAlik) 



(78) 



(79) 



and Eq. (I66p . we obtain for the scattering term /g'^''^ that governs the evolution of the statistically 
averaged quark distribution function the following expression: 

r(Ah\, \ 9 d f d'^k d^k' 



2N,dp~ J (2^)4 (2^)4 "^P(-^(fe - k')x){-ip'^k^g^ + ik,pn{SA';^{k)f;^^^{k',p)).{80) 

It is easy to see that Im[Ig] = because MJi(A;) = 6A*^"-{-k) and 5/^(^)(A;,p) = '5/*(^)(-A;, p). Let 
us calculate {5A1^{k)6f;^^^^{k' ,p)). Taking into account Eqs. (^-(^ and dUD we get 



{6A;{k,)6f;^^^{k2,p)) = (Mf °")(A.i)<5/;;«('="")(fc2,p)) 



+ 



(M^(^-*)(A:i)d/;;f ^"'■*)(A:2,p)) + (M;1(^'-*)( 



*a(s) , 



(81) 



where 



(M»(-")(^i)5/;f (^°")(fc2,p)) = 5 



Po-k2 — iO 



(M;i(^-*)(fci)<5/;j(;)(fc2,p)) = -9 [ d^p':^(e~')^.{.)(A;i)(<5/,"(^)(A;i,p')5/r^^^(fc2,p)), (83) 



2K 
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2E" 



Per ^2 ~ *0 



^^g^!'^^ (^"'^>Aw(^2)((^/°^^^(^ (84) 



Here we take into account that the correlator of the independent fluctuations is equal to zero. 

Now the calculation of the collision term Iq^^\x, p) is reduced to finding the explicit expressions 
for correlators in Eqs. (f82l) - (f84l) . Let us start with {5A"'^^°!^^\k)5A*^'^^°!^^\k')). One can easily get 

from Eq. (j75p (by multiplying the corresponding equation on 5A*°i^^^^\k') and taking the statis- 

tical average) that i5A°'}l^°!''\k)5A*"i^^°!'^\k')) is proportional to 8i\e^\{k)\) where |ef'^N(A;)| denotes 

the determinant of e'f\(k). Let us assume that is slowly varying in space-time 

function. Then the correlation function of collective classical field fluctuations is 

{5Al%"''\k)6A*;ff\k')) = {2^)H^'){k - A:05(lC)(fc)l)V(x)(^'^)- (^5) 



Here 



I^^i{k) flxes the strength of initial correlations of collective field modes and k^ = a;(^)(k) is the 



turbulent self-energy solution of the equation |e|'^-)(^)| = 0. 

The other correlators in Eqs. ([82]) - ([8l]) can be expressed through {Sfqj'^\k,p)6fq^^^\k',p')), 
where qi, qj are q,q or g. Because the spontaneous particle phase density fiuctuations are the 
solution of Eq. (|47|) . the corresponding correlation function should be proportional to 5{p^k^). 
The complete lack of initial correlations in velocities and positions means that the correlation 
function should contain 5*^'^)(p — p')6^'^\k — k'). Then, assuming statistical Poisson fluctuations 
(and, so, neglecting quantum statistics, etc.) and assuming, thereby, that the average value of the 
fluctuation in the squared number of particles in a certain volume V, {6N'^), is equal to the average 
number (A^) of particles in the volume V, we get 

{dfl^'\k,p)6f*f-'\k',p')) = i27rfi27ry%,EMx,p)6^^Hk - k')6^^\p - p')S{p>^k,). (87) 

Inserting expressions for correlation functions (|85|) and (|87p into Eqs. (|82|) - (j84|) and then into 
Eq. (j8ip . one can see that the collision term, Iq'^^\ in Eq. (j80p has the Fokker-Planck form with 
transport coefficients of drag, ^^(^.^(p), and diffusion, -D,y^(a;)(p): 

4^'^(->p) = ^ (^.Mxo(p)^^J^) + 1; Km(p)/.(->p)) • (88) 

Here 

^..w(p)=i^SS(P)+i^?;(5(p), (89) 

^t'™(x')(p) quasilinear diffusion coefficients that appear due to collective field fiuctuations, 

and -D^^"^*^ (p) are the diffusion coefficients due to particle fluctuations. It is noteworthy that in the 
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quasilinear approximation the collective field fluctuations contribute only to the diffusion coefficients 
of the Fokker-Planck equation ^28]. The quasilinear diffusion terms describe the isotropization 
processes in momentum space, while the generalized Balesku-Lenard terms [35] D^^"^*^p), A,y(^^-^{p) 
lead to thermalization and (local) equilibrium - the next stage after isotropizationlj In general, 
the Balescu-Lenard thermalization time is larger than the time scale of isotropization. Then, for 
relevant time scales, 

r)(..p) . A («i"'2,(P)^) . (90) 

where the quasilinear diffusion coefficients in the Abelian approximation are 

= I + (-<-■ + *(IC,(m)'„',.,(M). (91) 

4 Concluding remarks 

The isotropization time is the typical time at which the quark-gluon distributions become locally 
isotropic in momentum space due to diffusion in momentum space, and an important question is 
under what conditions isotropization of wQGP can be reached within time scales of relativistic 
heavy ion collisions. It is noteworthy that the value of isotropization tim^ depends not only on 
the instability rate but also on the initial value of the diffusion coefficient, and, thereby, on the 
initial value of the correlation function of collective field fluctuations. 

How the diffusion leads to the spreading out of the width of the distribution can be seen from 
the following simple mathematical example. Let us consider the diffusion equation 

(92) 



The solution can be written as 



where the Green's function 



dt duj"^ 



duj'G{Lo,Lo',t)f{io',to), (93) 



represents the function that obeys Eq. ([92]) and equals 5(a; — to') at t = to, and 

t 

d{t) = [ dt'D{t'). (95) 



to 



®It is noteworthy that the kinetic equation with Balescu-Lenard colhsion terms (as well as with Landau or Boltz- 
mann ones) is, in the thermodynamic sense, a kinetic equation of ideal gas: interaction contributes to dissipative 
quantities and does not contribute to thermodynamic ones. Then the Balescu-Lenard collision terms are, streakily 
speaking, improper near the equilibrium state if the state is far from the ideal gas one. 

^Because the corresponding collision terms depend, in general, on the nonAbelian dynamics of the gauge fields 
even at weak coupling, the value of isotropization time can be calculated only by means of numerical methods. 
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Then the width of f{io,t) is determined by the time dependence of the diffusion coefficient D{t). 
If, e.g., D{t) = D(io)exp(7(i — to)), then d{t) = (exp(7(t — to)) — 1) and for 7 > increases 
with time. 

Then, keeping in mind relations ([36l) . ([38]) . ([85]) . and (f9T]) between the diffusion term and the 
correlator of the fluctuations, we can suppose that, if energy stored in the initial collective fields' 
fluctuations is small (say, -D(to) is small), then even for strong initial anisotropy (high values of 7) 
during some initial transient time the width of the distribution function will change rather slowly 
with time (d(t) is small), and only after growth of collective field modes will the momentum width 
of distribution function start to grow. On the other hand, if the initial collective fields' fluctuations 
are strong enough, then the momentum width of the distribution function can grow rapidly from 
the very beginning, and in this case the system can reach the isotropical state while there are no 
developed instabilities with rapid growth of the fields, as have been noted in Ref. [20j . 

Note that any initial fiuctuations of the color fields are determined by fluctuations of the sources, 
and in the CGC approach fiuctuations of the sources result from fiuctuations in the color charge 
density in each of the colliding nuclei. Also, violations of boost invariance (really nucleus is not 
contracted into an infinitely thin sheet as was assumed in the original McLerran-Venugopalan 
model ^36j) result in the strengthening of fiuctuations of the color fields' sources. As a result of 
these perturbations extremely disordered color field configurations [37] can appear. Then, if strong 
fiuctuations take place, they, according to our analysis, generate effective collision terms that would 
result in fast local isotropization in momentum space of statistically averaged particle phase-space 
densitiesjfl A careful analysis of the initial conditions after a heavy ion collision would thus be a 
key ingredient in understanding the process of isotropization. 
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